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ABSTRACT 

We examine the prospects for studying the pre-reionization intergalactic medium (IGM) 
through the so-called 21 cm forest in spectra of bright high-redshift radio sources. We first 
compute the evolution of the mean optical depth r for models that include X-ray heating of the 
IGM gas, Wouthuysen-Field coupling, and reionization. Under most circumstances, the spin 
temperature Ts grows large well before reionization begins in earnest. As a result, r < 10^'^ 
throughout most of reionization, and background sources must sit well beyond the reioniza- 
tion surface in order to experience measurable absorption. HII regions produce relatively large 
"transmission gaps" and may therefore still be observable during the early stages of reioniza- 
tion. Absorption from sheets and filaments in the cosmic web fades once Ts becomes large 
and should be rare during reionization. Minihalos can produce strong (albeit narrow) absorp- 
tion features. Measuring their abundance would yield useful limits on the strength of feedback 
processes in the IGM as well as their effect on reionization. 

Key words: cosmology: theory - intergalactic medium - radio lines 



1 INTRODUCTION 



The "twilight zone" of structure formation, when the first galax- 
ies and quasars formed and eventually reionized the Universe, is 
the last frontier of galaxy formation studies. Observations are now 
beginning to probe this epo ch, through both direc t observations of 
high- redshift galaxies fe.g.. lBouwens et alj|2004t 
I2OO5I) and indirect constraints on reionization (e.g.. 



Taniguch^^U 
Fan et allZOOlf 



>• ■ Malhotra & Rhoadsll^ol iFan et alj|2005t IPage et a l 200g). But 



3P 
OOff). 



the chief lesson of these initial forays has been how little we know 
about the first luminous sources, their formation processes, their 
■ feedback on t he intergalac tic mediu m (IGM), and reionization it- 
' self (see, e.g.. lBarkana &L oeb 200ll: ICiardi'& Ferrara 2005). 

Thus a great deal of attention is now focused on develop- 
ing new approaches to observing the 2: > 10 Universe. Perhaps 
the most exciting is the 21 cm transition of neutral hydrogen (see 
iFurlanetto et all2006bl and references therein). Of the many appli- 
cations of this transition, the most popular possibility is to perform 
"21 cm tomography" of the IGM by mapping the redshift-space 
fluctuations in its brightness temper ature ( JIogan_&.Rees. J^7£; 
IScott & Reesll99Cl:lMadau et alll997l:IZaldarriaga et alJ2004 . Bv 
measuring the entire history of reionization, as well as its spatial 
structure, such observations could provide extraordinarily power- 
ful constraints on reionization and early structure formation. Un- 
fortunately, this technique also suffers from three intrinsic difficul- 
ties. First, it requires the spin temperature of the IGM, Ts, to dif- 
fer from the cosmic microwave background (CMB) temperature 
T.y. Otherwise, the two blend together and cannot be differenti- 
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ated. Unfortunately, this condition is not necessarily satisfied in the 
earliest phases of structure formation. Second, given realistic tele- 
scope parameters, tomography cannot resolve structures smaller 
than ~ 1 Mpc, which is well above the size of collapsed (or col- 
lapsing) structures ai z > 10. It therefore does not probe structure 
formation across all scales. Finally, tomography presents a number 
of truly formidable observational challenges, such as foreground 
contamination, ionospheric distortions, and terrestrial interference. 

In this paper, we will consider a complementar y probe im- 
mune to most of these problems: the "21 cm forest" JCarilli et alj 
2002; Furlanetto & Loeb 2002). The technique is the exact analog 
of the Lya forest that is so useful for studying the IGM at 2; < 6. 
Neutral hydrogen along the line of sight to a distant radio-loud 
quasar resonantly absorbs continuum radiation that redshifts into 
the local 21 cm transition. Thus a high resolution, high signal-to- 
noise spectrum can map the fluctuating neutral hydrogen field along 
this line of sight. Interestingly, the 21 cm forest has a surprisingly 
long historv: .Field Ll959a') attempted to observe it toward Cygnus 
A several years before the Lya forest was even discovered (in the 
process, he ruled out an fi = 1 universe full of neutral hydrogen). 

The 21 cm forest bypasses all of the difficulties with tomo- 
graphic measurements. First, the CMB no longer acts as the back- 
ground source, so the IGM remains visible even if Ts = T^. 
Second, the IGM can be resolved spatially on fine scales (tens 
of kpc) given a sufficiently bright background source, allowing us 
to probe individual filaments and minihalos with square-kilometer 
class telescopes. Finally, spectroscopy of bright point sources is rel- 
atively straightforward from an observer's perspective, and 21 cm 
forest observations are immune to many of the systematics that be- 
set tomography. 

We expect to find four types of absorption features. First, the 
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mean level of absor ption measures the global evolution of the ra- 
diation background iCarilli et alJl2002l) . Second, sheets and fila- 
ments in the "proto-cosmic web" produce stronger absorption fea- 
tures, just as they do in the Lya forest: these features provide in - 
sight into the first phases of structure formation iCarilli et alj2002l) . 
Third, a fair fraction of the IGM gas may condense into "miniha- 
los" before reionization. These are virialized neutral clumps that 
are too small to fragment a nd for m stars bu t may have important im- 
plications for reio nization iOh & Haimanl l2003l : llliev et alj|2005d : 
[Ciardi et al. 2 0051) . No other method has been devised to detect 
them directly, but they pr oduce relatively strong 21 cm absorption 
JFurlanetto & L oeb 2002). Finally, HII regions produce transmis- 
sion gaps in the absorption spectrum; their evolution would obvi- 
ously constrain the evolution of the ionizing background as well as 
its spatial structure. While the first three of these classes have been 
considered before, they have not been examined in the context of 
self-consistent histories of the 21 cm background. The purpose of 
this paper is to re-examine each kind of absorber in light of this 
broader context. We will examine the features each imprints on the 
21 cm forest, the epochs at which they are important, and their ob- 
servable implications. 

The remainder of this paper is organized as follows. We de- 
scribe some basic properties of the 21 cm transition in iJ2| We then 
discuss each of the four types of spectral features in turn: the mean 
level in Sj3| HII regions in !J4| the cosmic web in [Js] and minihalos 
in |6| We conclude and discuss potential observations in ^ 

In our numerical calculations, we assume a cosmology with 
^m = 0.26, Qa = 0.74, 0,1, = 0.044, H = lOO/i km s"^ Mpc'^ 
(with h — 0.74), n = 0. 95, and erg = 0.8, consistent with the 
most recent measurements fSperael et al. 2006), although we have 
increased erg from the best-fit WMAP value to improve agreement 
with weak lensing. We quote all distances in comoving units, unless 
otherwise specified. 



uses the expansion rate to determine the total column. For a single 
isolated cloud, the central optical depth is 

1/2 

Tel ^ 0.076 ( ^^^ I ( ^^^^ _,;;"^_, I , (3) 



Tk 



Ts 1021 cm" 



where Tk is the kinetic temperature, A'hi is the total column den- 
sity of neutral hydrogen, and we have assumed pure thermal broad- 
ening. Thus, clouds must have column densities close to those of 
damped-Lyof absorbers in order to cause substantial 21 cm absorp- 
tion. 

The spin temperature Ts is determined by a competition be- 
tween three processes: atomic collision s, scattering of CMB pho- 
tons, and scattering of Lya photons iWouthuvsenI II952 : iField 
1958). In equilibrium. 



Tr^ = 



_/~, ~t~XQ<_/^ -t~jjc-t7> 



(4) 



1 -t" Xa I Xc 

Here Xc is the sum of the coUisional coupling coefficients for H- 
H interactions (Zvgelman 2005) and H-e^ collisions (Smith 196g; 
Liszt 2001); it is of course proportional to the local density. The 
middle term describes the Wouthuysen-Field effect, in which ab- 
sorption and re-emission of Lya photons mix the hyperfine states. 
The coupling coefficient is 



1.81 X 10"(l + z)"^S'„Ja, 



(5) 



where Sa is a factor of order unity describing the detailed atomic 
physics of t he scattering process ( Chen & Miralda-Escude 2004j 
lHiratal2005l) and Ja is the background flux at the Lya frequency in 
units cm~2 s~i Hz"'^ sr^^; the Wouthuysen-Field effect becomes 
efficient when there are ~ 0.1-1 photons per baryon near this fre- 
quency. It couples Ts to an effective color temperature Tc, typically 
fc ^ Tk (Field 1959b.). We use the numerical fits of .Hiratij^ 1,2001) 
for Sa and Tc. 



2 THE 21 CM TRANSITION 

We review the relevant characteristics of the 21 cm transition here; 
we refer the interested reader to Furlanetto et al. ( 2006b) for a more 
comprehensive discussion. The 21 cm optical depth of a patch of 
the IGM is LField.l959a.) 



l+z\3/2 



9.6 X W-^XHi{l + 5) (^^) 



Ts 



H{z)/{1 + 
d^li/drii 



(1) 



where xyii is the neutral fraction, 5 is the fractional overdensity, 
T~i is the CMB temperature at redshift z, Ts is the spin temper- 
ature, and du||/dr|| is the line of sight velocity in physical units. 
The brightness temperature of the patch viewed against a source 
of brightness temperature Trg is (neglecting peculiar velocities for 
simplicity) 



sn 



27a;Hi (1 + 5) 



Ts - r„ 



o.o23y V 



0.15 1 + 2X1/2 



Ts 



0.023/ \Q,rr,h? 10 



mK. 



(2) 



Note that &Tb oc T^g/Ts for Trg > Ts, the regime relevant for 
this paper. Thus the IGM temperature is crucial to estimating the 
signal strength. 

Equation Q applies only to the expanding IGM, because it 



3 THE AVERAGE SIGNAL 

We will first consider the evolution of the mean optical depth with 
time, which is determined purely by the mean neutral fraction 
xhi{z) and Ts{z); note that we neglect fluctuations in these quan- 
tities here. These depend on the thermal history of the IGM (includ- 
ing adiabatic expansion. X-ray heating, and shocks), the ultraviolet 
(UV) background (which determines the Lya coupling and hence 
spin temperature), and the ionizing efficiency. Unfortunately, none 
of these factors are well-constrained at high redshifts. Rather than 
attempt to make robust predictions, we will therefore content our- 
selves with examining a range of simple histories to develop some 
intuition about the qualitative features^ 

We will use the simple model o fJFurlanettd llOOd) to compute 
the optical depth history; we refer the reader there for a detailed 
discussion of the various prescriptions. In brief, this model assumes 
that stars dominate the radiation background and sets the total star 
formation rate via the rate at which gas collapses onto galaxies, 
d/coii/df, where /coii is the collapse fraction (or the mass fraction 
in halos with virial temperatures Tvir > 10* K). The ionizing effi- 
ciency is (" = Aiicf*fcscNion, where /* is the star formation effi- 
ciency, /eac is the escape fraction of ionizing photons, A^'ion is the 
number of ionizing photons produced per baryon incorporated into 
stars, and ^hc is a normalization constant accounting for helium in 
the IGM. 

The rate of X-ray heating is also assumed to be proportional 
to d/coii/di, with a fiducial proportionality constant fixed by the 
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Figure 1. Optical depth evolution in one-population reionization histories, (a): Spin temperature. The dotted curve shows T-y for reference, (b): Ionized fraction 
(c): Optical depth in neutral regions (d): Mean optical depth. The solid and dot-dashed curves use our fiducial Pop II model with fx = 1 and 0.2, respectively. 
The long- and short-dashed curves use the Pop III model with fx = 1 and 5, respectively. 



local relation between star formation r ate and X-ray lumin osity 
JOrimm et alJl2003l: iRanalli et alJl2003l:lGilfanov et 31112004 . We 
let fx measure the efficiency at high redshifts relative to this lo- 
cal value (measured between 0.2 a nd 10 keV). We include a sim- 
ple model for shock heating from iFurlanetto & Loea i2004) that 
describes shocks in the high-redshift cosmic web reasonably well 
jKuhlen et alJ20 0(^. 

The final element is the UV background responsible for the 
Wouthuysen-Field effect. We again assume that the emissivity is 
proportional to d/coii/d£ and use the (low-metallicity) Popula- 
tion II and (very massive ) Population II I (henceforth Pop II and 
Pop III) spectral fits of Barkana & Loeb (2005) to estimate the 
efficiency with which these photons are produced. Note that we 
properly incorporate higher L yman-series transitions fairatal2005l : 
IPritchard & Furlanettol2005h . 

Although there are obviously a number of free parameters in 
this model, the general features are easy to understand. The 21 cm 
background contains three major transitions: the points when Tk 
first exceeds T-y, when Lya coupling becomes efficient, and when 
reionization occurs. Because it affects all the backgrounds equally, 
/* simply shifts everything forward and backward in time. The X- 
ray efficiency fx only affects the kinetic temperature, while /esc 
only affects xhi. The stellar initial mass function - especially the 
choice between Pop II and very massive Pop III stellar populations 
- affects the ratio of Lya photons to ionizing photons and moves 
the onset of Lya coupling relative to reionization. Because they are 
so much hotter, Pop III stars push this transition forward in time, 
compressing the features in the 21 cm background. 

Figure Q shows some results for single-population models 
(without feedback). For Pop II stars, we take /* = 0.1, /esc = 0.1, 
Mon = 4000, and /x = 1 (solid curves). For Pop III stars, we 
take /* = 0.01, /esc = 0.1, fx = 1, Men = 30, 000, and fx = 1 
(long-dashed curves), for a somewhat smaller overall ionizing ef- 
ficiency. We also show results for Pop II stars with a low X-ray 



heating efficiency fx =0.2 (dot-dashed curves) and Pop III stars 
with a high X-ray efficiency fx~5 (short-dashed curves). 

Panels a and b show Ts{z) and Xi(z) = 1 — x-h_i{z). With 
this normalization, reionization ends at 2 ~ 7. The most important 
point is that Ts increases beyond T-y well before reionization com- 
pletes, although the transition does occur significantly later for Pop 
III stars. Crucially, Ts also saturates at a smaller value for these 
stars. 

Panel c shows the average optical depth in neutral regions (i.e., 
the contrast between neutral absorbing regions and empty HII re- 
gions). Panel d shows the mean optical depth as a function of red- 
shift (i.e., the optical depth averaged over long stretches of the spec- 
trum, including both neutral and ionized regions). At sufficiently 
high redshifts, Ts « T^ and f evolves slowly. Once Lya cou- 
pling becomes strong, f increases somewhat because (in all of these 
models) X-ray heating kicks in somewhat later. However, it does 
not take long for X-ray heating to dominate the temperature. From 
this point, the optical depth decreases rapidly: for our fiducial Pop 
II parameters, r < 10~^ almost before reionization begins. Even 
though massive Pop III stars produce comparable kinetic temper- 
atures, the IGM optical depth remains about an order magnitude 
larger if they dominate. This is purely a function of the spin temper- 
ature (and lack of Lya photons), so it would also occur if quasars 
dominated the radiation background. Thus the 21 cm forest may be 
a more powerful probe if Pop III stars, or some other hard source, 
dominate. Note that emission against the CMB does not depend 
on Ts so long as Tg ^ T-,. Thus comparing the mean levels of 
emission and absorption could help to identify the type of sources 
during reionization. 

In any case, FigureQillustrates a crucial point about the 21 cm 
forest: observing high optical depth features requires sources that 
shine well before reionization. At 2 > 15, r can indeed be quite 
large - nearly 10% in some models - but during reionization it is 
likely to be up to three orders of magnitude smaller. Unfortunately, 
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Figure 2. As FigureQ but for models that include photoheating feedback. The solid and dot-dashed curves use our fiducial Pop II parameters. The long and 
short-dashed curves assume that Pop III (Pop II) stars form in cold (hot) regions. In each case, the two curves take different feedback speeds (see text). Note 
that the "double reionization" seen in the long-dashed curve is difficult to achieve in realistic situations. 



this presents a quandary: radio sources are likely to be associated 
with ionizing sources, so they too will be extremely rare at such 
high redshifts. The 21 cm forest has the inherent difficult that the 
strongest absorption can only occur for the most distant sources, 
while it may be nearly unobservable for the most common sources 
during reionization. 

Of course, feedback (among other processes) could easily 
cause more complicated histories. Figure|2|shows some examples. 
The solid and dot-dashed curves assume our Pop 11 parameters but 
include strong photoheating feedba ck^ (with varying speeds; see 
[Furlanetto & Loeb 2005 and Furlan ettci2006l for details). The long 
and short-dashed curves increase the effects of feedback by assum- 
ing that photoheating is accompanied by a simultaneous change in 
stellar population; again we allow two different feedback speeds. 
The parameters for each population take our fiducial values, ex- 
cept that we set /* = 0.1 for Pop 111 stars in order to exaggerate 
the decline in ionizing efficiency. Note that the long-dashed curve 
yields a short period in wh ich recombinations domi nate. This is 
unlikely to occur in reality JFurlanetto & Loebll2005l) . as we have 
stretched several parameters to achieve it (raised the Jeans mass in 
photoheated region to its maximal value, assumed unrealistically 
fast feedback, and conflated the transition between stellar popula- 
tions with photoheating). But it helps to illustrate how such an era 
might affect the signal. 

The principal effect of feedback is to prolong reionization by 
Aa ~ 5-10. One might then hope for a long period of strong con- 
trast between neutral and ionized regions. But panel c shows that 
the optical depth in neutral gas nevertheless falls to ~ 10""' by 



^ Specifically, we assume that photoheating increases the minimum galaxy 
mass to Tvir = 2 X 10^ K (Thoul & Weinbera 1996; Kitavama & Ikeuchi 
I2OO0I). even though recent work suggests photoheating is actually less effi- 
cient iDiiksfi-a et all2004l) . 



z ~ 10 in all cases. This is because X-ray heating, and Lya cou- 
pling, establish themselves in the early stages of reionization, when 
feedback is irrelevant. Thus for a fixed endpoint to reionization, 
feedback weakens the 21 cm forest by pushing the X-ray and Lya 
coupling transitions farther back in time. Again, the prospects are 
somewhat better if Pop 111 stars dominate the early stages. 

Now that we have some idea of the optical depths involved, it 
is useful to estimate their observability. The minimum background 
source flux density required to detect an absorption feature with 
signal to noise S/N is 



16 mJy 



/ S/N 10~^ 10*^ m^ %, 



^5 r yloff 400 K 
/ 100 kHz lweek\^/^ 



V Ai/ch iint 



-} 



(6) 



where A^s is the effective area of the telescope, T^ys is the sys- 
tem temperature, Afch is the bandwidth of each channel (assumed 
to be narrower than the feature of interest), and tint is the total 
"on-source" integration time. Here we have assumed simultaneous 
observation of two orthogonal polarizations, which is normal for 
dipole antennae, and taken telescope parameters sim ilar to those 
expec ted for the SKA observing at 2 = 10 (see also lCarilli et alJ 
I2OO2I) . At these low frequencies, Tsys is dominated by the Galac- 
tic synchrotron background; on a quiet portion of the sky, this is 
typically (e.g.. lFurlanetto et alj2006bh 



'sky 



180 f '^ y 

V 180 MHz/ 



K. 



(7) 



As we shall see, observing any of the absorbers we will discuss re- 
quires rather luminous sources. For example, the 2 = 10 equivalent 
of Cygnus A would have 5" < 20 mJy 1 Carilli et al. 2002) . 
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Figure 3. Differential abundance of HIl regions intersected per unit redshift 
and per unit radial width at z = 10. The dotted, dot-dashed, short-dashed, 
long-dashed, and solid lines take Xi = 0.1, 0.3, 0.5, 0.7, and 0.9, respec- 
tively. 



Figure 4. Differential abundance of Hll regions intersected per unit redshift 
and per unit radial width at 2 = 15. The solid, dotted, and dashed lines take 

Xi = 0.05, 0.1, and 0.5, respectively. The thick and thin curves assume 

2/3 



C, ex m^ and m^ 



4 HII REGIONS 

Ionized bubbles create gaps in the absorption. In this section, 
we will estimate the abundance and sizes of such transmission 
gaps throughout reionization. Carilli et al. 1^02) implicitly in- 
cluded HII regions in their simulation of the 21 cm forest. How- 
ever, their box was only Ah ~^ Mpc across. Recent analyti c mod- 
els i Furlanetto et alj l2004h and simulations I Iliev et al] l2005iJ: 
IZahn et al . 2006) have shown that ionized bubbles rapidly exceed 
this size, so much larger boxes - or analytic models ~ are required 
to sample the bubble distribution reliably. 

We will therefore use a simple analytic model for the bub- 
ble size distribution I Furlanetto et al. 2004, 2006a) that provides 
a good match to simulations iZahnetal. 2006). The model as- 
sumes that the distribution of ionized material is driven primarily 
by the strong clustering of ionizing sources. We will ignore recom- 
binations (or more precisely their spatial fluctuations, which ap- 
pears to be a good assumption during the early and middle stages 
of reionization for which the 21 cm absorption signal is large; 
IFurlanetto & Oh-2005.1 . In that case, the mass that any galaxy can 
ionize is simply ^mgai. Thus, an isolated region is fully ionized if 
/coll > 1/C; because the collapse fraction is itself a function of the 
large-scale den sity field, this con dition can be used to specify the 
ionization field. IFurlanetto et alJi2004) showed that the comoving 
number density of ionized bubbles ni,{m) is then 



mnt{m) = \ 



1£ 

■K m 



diner 



dlnm 



Bo 

a{m) 



exp 



B'^{m,z) 
' 2cr2(m) 



(8) 



where Bo and B are density thresholds from the condition on 
the collapse fraction, a^{m) is the variance of the density field 
smoothed on mass m, and p is the mean comoving matter density. 
Equation jsj can easily be gener alized to allow the ioniz ing effi- 
ciency to vary with galaxy mass iFurlanetto et alJ l2006a ). so that 
C = C{mh) and irih is the halo mass. As we increase the im- 



portance of massive galaxies, the bubbles grow larger, especially 
during the early phases of reionization. 

With equation (Sj, the number density of transmission gaps in 
a 21 cm forest spectrum is 



dn(> Ai/obs) dr / 2 / . 
-j ^ d~ I <imn6(m)7rrcs(m, Ai^obs 



(9) 



where r-cs is the maximum impact parameter for a bubble of 
mass m to have a spectral length (in observer's units) of Ai^obs 
and dr/dz is the comoving line element. For reference, Ar ~ 
1.7(A!/obs/100 kHz) Mpc at 2 = 10. 

Figure|3|shows the resulting distribution (transformed into dif- 
ferential form) for a variety of ionized fractions at z = 10 (we have 
adjusted (" in each case so as to recover the desired Xi). Except at 
the smallest Xi, the distribution is fairly flat until it cuts off sharply 
at relatively large widths. This is because ni,{m) has a well-defined 
characteristic size Rf,, corresponding to the cutoff. Narrower trans- 
mission features arise from a mixture of smaller bubbles and lines 
of sight that pass near the edges of bubbles with R ~ Ri,. Be- 
cause Rt rapidly grows from a few tenths to a few tens of Mpc 
iFurlanetto et aL.2004.) . the transmission gaps are relatively large 
spectral features. This contrasts with a model in which galaxy HII 
regions remain isolated, which would have gaps at most a few tens 
of kHz wide: thus galaxy clustering is extremely important for the 
distribution of these bubble features (and hence the 21 cm forest 
provides a sharp test). Another key feature is that the total num- 
ber density of distinct gaps remains roughly constant with Xi (at 
least until Xi > 0.75). This is because most new ionizations oc- 
cur through the growth of existing bubbles (which also merge with 
their neighbors) rather than through the formation of new bubbles 
( iFurlanetto & Oh 2005). As a result, the most obvious bubble fluc- 
tuations will actually occur relatively early in reionization, when 
bubbles do not carpet the spectrum. 

Of course, the amplitude of the transmission feature is de- 
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termined by the absorbing neutral gas surrounding tlie HII re- 
gion (e.g., FigureQ;). For our fiducial Pop II stars, it drops from 
T ~ 10"^ to ~ 10"* over the range Xi ~ 0.1-0.4. Clearly, the 
contrast will be small even in the early phases of reionization. Thus, 
by equation J6}, observing HII regions would require a source of 
comparable luminosity to Cygnus A. On the other hand, if massive 
Pop III stars dominate, the mean optical depth could be an order 
of magnitude larger (see FigureQ:). In this case, Mpc-scale HII re- 
gions could be observed against mjy radio sources, which may be 
reasonably common at z < 12 (see iQ. 

Figure|4|shows how the size distribution changes with the pa- 
rameters of the ionizing sources. The thick curves show an iden- 
tical model to Figure |3| except at z — 15. The thin curves as- 
sume that ^ P C m/ : thus massive ga laxies shift the HII regions to 
larger scales ( Furlanetto e t al. 2006a). This manifests itself in the 
spectrum by moving the radial cutoff to slightly larger frequency 
widths; the effect is largest at small Xi, where the bubbles will be 
easiest to observe. Obviously the range of dn/dz is relatively small 
in these models, but measurements will still help to constrain more 
complicated models (including feedback, for example) to test the 
basic premises of the Furlanetto et al. 1 2004) model, and to mea- 
sure Xi{z). 

Finally, a comparison of the Xi =0.1 and 0.5 curves in Fig- 
ures|3|and|4|shows that the size distribution of HII regions at a given 
ionized fraction does not depend strongly on redshift. The shapes 
are nearly identical; the normalization changes primarily through 
the line element dr/dz. 



thermal broadening (additional sources can only decrease the peak 
r) and Ts ~ Tk (valid at high densities and/or for a substan- 
tial soft-UV background), we find that the mean overdensity corre- 
sponding to a feature of optical depth r is 



1 + S 



94 



f^V f^^v (j^) 

Vo.oiy Viook; \i + zJ 



0.15 



0.023 



(11) 



Overdensities ~ 100 only occur in virialized objects. FromFigs.^ 
and|2lj, we see that in most cases Ts > 100 K long before reioniza- 
tion enters its rapid phase, even if X-ray heating is relatively weak. 
Clearly, in these sorts of models cosmic web absorption features 
will be extremely rare during reionization. Feedback only makes 
the prospects worse by pushing heating even further back in time 
relative to the end of reionization. 

Recall that Carilli et al. ( 2003) found significant absorption 
even at z ~ 10, which may appear to disagree with our claims. 
This is because their simulation ignored X-ray heating so had 
Ts ~ 30 K throughout much of the IGM at 2: = 10. Our models 
show that X-rays can easily surpass this threshold by z ~ 12-15. 
On the other hand, they had Ts > 100 K everywhere by z = 8 - 
and indeed the number density of absorbers fell by a factor of ten 
over that time interval. So our principal conclusion, that Ts must 
remain small in order for sheets and filaments to be visible, agrees 
with their simulation. This just illustrates the crucial importance of 
considering the signal in the broadest possible context. 



5 THE COSMIC WEB 

A third type of feature is the exact analog of the Lya forest: the 
fluctuating absorption from sheets and filaments beginning to con- 
dense out of the IGM. Thes e can appear as nar row r > 1% absorp- 
tion spikes in the spectrum (Carilli et al. 2002). Naively, one might 
expect that these features would strengthen with time as the cos- 
mic web accumulates more mass through ongoing structure forma- 
tion. However, as usual we must also consider the broader context 
of the radiation backgrounds. These increase Ts everywhere and 
hence decrease the strength of absorption even inside s heets and fil- 
ament s. For example, in the simulations examined bv lCarilli et alj 
y003), the number density of absorbers with r > 0.02 actually 
decreased from ~ 50 to ~ 4 over the range 2 = 10-8. 

To study this problem in a more general context, we follow 
ISchavd |200l|) and consider a cloud with characteristic hydrogen 
density uh that is in (or near) hydrostatic equilibrium. In such 
a cloud, dP/dr = —G p M{r)/r^, where P ~ cj p is the gas 
pressure, Cs is the sound speed, and Ad{r) is the mass enclosed 
within a radius r. Hydrostatic equilibrium thus sets a characteris- 
tic "Jeans length" Lj for the cloud through c^ p/Lj ^ G p^ Lj, 
which describes the typical comoving length scale over which the 
cloud achieves a given overdensity. 



Lj 



7.1(1+5) 



-1/2 



10 0.15 



Tk 

100 K 1 + 2 n™^2 



1/2 



kpc.(lO) 



Despite its obvious approximations and shortcomings, this simple 
model does a surprisingly good job describing the Lya forest at 
2 ~ 3, which is composed o f the same sheets and filaments as the 
21 cm forest iSchavell200ll) . We therefore expect it to provide a 
reasonable description of high-redshift features as well. 

The column density of the cloud is then A^hi ^ nn Lj, and 
we can compute its optical depth via equation J3j- Assuming pure 



6 MINIHALO ABSORPTION 

A final set of features come from minihalos. These gas clouds can 
have extremely large column densit ies and hence produce rela- 
tively strong absorption. Furlanetto & Loebl J2002) computed the 
absorption profiles of minihalos and their abundance for a particu- 
lar thermal history (based on the simulation of C arilli et al. 2003) . 
As above, our goal here is simply to consider the signal in more 
general models. 

We assume that min ihalos form with a comovin g number den- 
sity nh( m) given by the Press & Sch echten il97 4) mass function 
(using the lSheth & Tormen 1999 mass function affects our results 
by < 25%). We assume that the gas settles into an isothermal dis- 
tribution (at the virial temperature) inside the dark matter potential, 
which follows the Navarro et al. 1 1997) distribution. Other gas pro- 
files will modify the absorption statistics, though not dramatically. 
The excess optical depth al ong a line of sight with im pact parame- 
ter a through a minihalo is JFurlanetto & Loebl2002l) 



T^,i,{v) = 1.91 X 10*^ 



dR 



Ts{r)b 



exp 



-v\u) 



b2 



(12) 



where h^ = 2fcsrvir/m.p is the Doppler parameter, v{u) — c(u — 
i^o)/i^2i, I' is the frequency at the minihalo (without the cosmolog- 
ical redshifting to the observer), 1/21 = 1420 MHz, r^ = a^ +R^, 



R„ 



, and R is measured in physical kiloparsecs. Here 



we have neglected the infall region around eac h minihalo, which 
can a dd a significant amount of absorption 1 Furlanetto & Loebl 
I2OO2I) but is difficult to model accurately. Note that, because the 
density varies with radius, Ts — Ts{r) even though the gas is 
isothermal. The peak optical depth can be large - tens of percent 
- though lines of sight with Q>0.1rvir usually have r < 0.02 
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JFurlanetto & Loebll2002h . The typical line widths (resulting from 
thermal broadening) are A/^obs ~ 2 kHz. 

As with filaments, the UV background will affect Ts inside 
minihalos. However, in this case it has only a modest effect because 
the local densities and temperatures are large, so that Ts ~ Tk 
is a reasonable approximation everywhere. Nevertheless, even if 
id 3> 1 everywhere, the minihalo optical depth greatly exceeds 
that of the IGM, because of their high column densities. In such 
a situation, the only difference is that equation <12> then measures 
the excess optical depth over the baseline IGM level (at least to 
zeroth order). This mark s an important contra st with attempts to 
study minihalo emission (Iliev et al. 2002, 2003): even when IGM 
emission overwhelms minihalos ( Fu rlanetto & Oh 2006') . they are 
still distinguishable in absorption through their large columns. 

The abundance of minihalo features is then 

dn(> r^h) dr /""'' , ^ 2. ^ ,,^. 
-T ^ 'a~ I amnh[m)Tvr^{m,Tn-,h), (13) 

where rT(Tn,rinh) is the maximum impact parameter (in comov- 
ing units) to have an optical depth greater than Tmh- The maxi- 
mum minihalo mass 77x4 is assumed to have Tvir = 10^ K, where 
atomic line cooling b ecomes efficient and stars can form (e.g., 
iBarkana & Loebll200ll) . The minimum mass is more problematic. 
The Jeans mass obvious ly increases as X-ray s (or other processes) 
heat the IGM (see, e.g.. lGnedin & Hui' 1998). But the quantitative 
effects of thermal feedback remain controversial. Oh & Haimaa 
i2003h argued that X-rays produce an "entropy floor" that effi- 
ciently suppresses accretion onto minihalos, and they showed that 
21 cm absorption becomes much weaker in such a scenario. How- 
ever, if the energy is injected after the gas has already achieved 
relatively high densities. X-ray heating may actually have only a 
modest effect on minihalo formation. Perhaps the most interest- 
ing implication of minihalo measurements will be to constrain this 
feedback, so we will consider a range of efficiencies. For the sake 
of transparency, we simply assign the minimum mass by the Jeans 
mass corresponding to a particular IGM temperature. 

Figure |5| shows the number density of minihalo absorbers at 
2: = 10 and 20 (solid and dashed lines, respectively). At each 
redshift, the top curves assume Tk = T^diz) (the temperature 
in the abse nce of structure fo rmation processes, computed from 
RECFASTi fSeager et alJll999^ . The remaining curves take Tk = 
20, 100, and 1000 K, from top to bottom within each set. In a cold 
IGM, the number of minihalo features is large (with > 100 ab- 
sorbers per unit redshift in this r range at 2 = 10) and they can 
be rather strong (with >20 r > 0.1 absorbers per unit redshift). 
There are about an order of magnitude fewer features at 2: = 20, 
simply because structure formation is so much less advanced at that 
time. 

But the number of absorbers decreases rapidly as Tk in- 
creases, especially for large optical depths. This is because the 
smallest minihalos produce the strongest features: by equation J3}, 
Tnih oc Ayirryir/Tg PC T~} , whcrc Ayjr is the mean overdensity 
in virialized objects iBarkana & Loebl200ir) . Thus raising the IGM 
temperature by even a small amount eliminates the strongest fea- 
tures. On the other hand, massive minihalos produce the more com- 
mon T ~ 0.01 absorption features, so they remain even if thermal 
feedback is strong. The sensitivity of the (strong) minihalo features 
tells us that 21 cm absorption can indeed be used to quantify the ef- 
ficiency of thermal feedback and thus to constrain models of early 
structure formation. If these objects continue to form throughout 
the reionization era, they may therefore present the best possibility 
for strong features in the 21 cm forest. 
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Figure 5. Differential number density of minihaio absorption features at 
z = 10 (solid curves) and 2 = 20 (dashed curves). From top to bottom, 
each set takes the minimum minihalo mass to be that of a uniform medium 
at Tk = TadCz), 20, 100 and 1000 K. Here T^dCz) is the temperature 
without including any heating from structure formation. Note that, at 2 = 
20, the Tk = 10^ K curve lies below the plot range. 



7 DISCUSSION 

We have calculated four kinds of "21 cm forest" absorption sig- 
natures. First, we computed the evolution of the mean optical 
depth in a set of plausible histories of early structure formation. 
We found that the optical depth can be relatively large (r > 0.03) 
at sufficiently early epochs, but it begins to decline rapidly once 
X-ray heating and Lya coupling both become efficient (because 
T (X Tg^). By the time reionization begins in earnest, r ~ 10^'' if 
Pop II stars dominate the radiation background. The optical depth 
can be somewhat higher if very massive Pop III stars dominate, 
because they cause weaker Lya coupling. 

Thus the best bet for observing the 21 cm forest is well be- 
fore reionization. HII regions may be the easiest features to ob- 
serve, because they are relatively large ( > 10-100 kHz) even in 
these early stages (although the contrast with neutral gas may still 
be relatively small). Their detection would allow us to measure 
the size distribution of ionized bubbles, a key quantity in under- 
standing reionization. Sheets and filaments in the cosmic web, as 
well as minihalos, produce strong but narrow ( ^ 2 kHz) features. 
The vast majority of cosmic web features most likely fade well 
before reionization, because only virial overdensities can produce 
strong absorption once Ts ^ 100 K. Minihalos may persist longer, 
because they have much higher column densities. We expect that up 
to ~ 100 absorbers with r > 0.01 may be visible per unit redshift 
at 2 = 10 (or ~ 8 at 2 = 20), and their abundance provi des a sensi- 
tive measure of thermal feedback iOh & Ha iman'20Q3>. They may 
also provide insight into the small -scale ma tter power spectrum, if 
the astrophysics can be separated Jlliev et a l. 2002, 2003). Regard- 
less, the 21 cm forest avoids the confusion between miniha los and 
the IG M that plagues tomography measurements 1 Furlanett o & OhI 
l200d) . 

We have considered each of these absorber classes in isola- 
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tion, but of course in reality tliey will be intermingled. Sheets, fila- 
ments, and minihalos will fluctuate around the baseline absorption 
provided by the mean signal, while HII regions will appear as gaps 
(although miniha los may still appear in such regions un til they are 
photoevaporated: IShapiro et all2004l : llliev et al]l2005bl) . This will 
actually make the ionized bubbles easier to detect, because the gaps 
will be more obvious even if the mean level of absorption is small 
(i.e., both the mean level and variance of absorption will decrease 
in HII regions). On the other hand, it may be somewhat difficult to 
separate minihalo and cosmic web absorbers, because both produce 
narrow, strong lines. 

The 21 cm forest is a much simpler observation than large- 
scale 21 cm tomography, because we can essentially ignore all of 
the wide-field imaging and foreground problems that plague those 
observations. Thus the feasibility of observing the 21 cm forest 
boils down to two issues: raw sensitivity and the existence of back- 
ground sources. Equation J6j shows that the required source bright- 
ness is rather large. As a result, large collecting areas - comparable 
to the Square Kilometer Array - are required for meaningful con- 
straints. Somewhat fainter objects can be used to detect fluctuations 
statistically (for example, minihalos and the cosmic web provide an 
extra source of apparent noise); however, such games will only buy 
a fac tor of few in source flux, so the limits are still relatively strin- 
gent JCarilli et all2002h . 

The most promising candidates for such bright objects 
(S > 6 mjy) are radio-loud quasars; gam ma-ray bursts and hyper - 
novae most likely do not reach mJy fluxes Hoka & Meszaroa2005l) . 
Do bright radio-loud quasars exist at high redshifts? Assuming that 
their number density traces that of optical qu asars, the sky coul d 
contain ~ 1000 (10) sources at z > 8 (12) jCarilli et al. 2002). 
More detailed models, calibrated to the lower-redshift radio lu- 
minosity function, predict ~ 2000 sources acros s the sky with 
S > Q mJy at 8 < z < 12 faaiman e t ai]|2004 . although this 
prediction is quite sensitive to the assumed radio-loudness parame- 
ter. Thus it appears quite reasonable to expect a fairly large number 
of appropriate lines of sight with 2; < 12. Intriguingly, such bright 
sources could alread y have been detecte d in higher frequency sur- 
veys, such as FIRST JBecker et alJl99j) . The difficulty lies in sep- 
arating t hem from the much more numerous low-redshift interlop- 
ers (e.g.. lHaiman et all2004) . 

Unfortunately, at extremely high redshifts (2 > 12) we will 
likely be limited to just a few lines of sight. This is particularly 
unfortunate because our models show that high optical-depth re- 
gions most likely only exist at such epochs. X-ray heating and 
Wouthuysen-Field coupling genera lly precede reionization by a 
long time interval iFurlanetta l200(^. decreasing the optical depth 
during reionization. Thus the rapidly declining source density and 
the increasing intrinsic signal work against each other, which may 
ultimately limit the usefulness of the 21 cm forest. Of course, 
this is not coincidental: bright quasars are most likely associated 
with massive galaxies, which in turn must be rare until reioniza- 
tion gets underway. Nevertheless, even a few high-redshift lines 
of sight would be true prizes, because they would allow us to ob- 
serve the first phases of structure formation with unprecedented 
spatial resolution and to constrain the radiation backgrounds in 
the early Universe. In this context, the recent downward revision 
of the WMAP optical depth (Spergel et al. 2006; Page et al. 2006) 
is certainly good news, because it suggests that heating and spin 
temperature coupling may have occurred at somewhat later times, 
and the requirements on background sources may be somewhat less 
stringent. Given the many uncertainties about the 2^6 Universe, 
lower-redshift sources will also offer meaningful constraints (for 



example, ruling out weak spin-temperature coupling from massive 
Pop III stars and testing feedback models for minihalos), even if 
the most plausible cases have considerably smaller signals. 

I thank the Tapir group at Caltech for hospitality while some 
of this work was completed. 



REFERENCES 

Barkana, R., & Loeb, A. 2001, Physics Reports , 349, 125 

— . 2005, ApJ , 626, 1 

Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ , 450, 559 

Bouwens, R. J., et al. 2004, ApJ ,616, L79 

Carilli, C. L., Gnedin, N. Y., & Owen, F. 2002, ApJ , 577, 22 

Chen, X., & Miralda-Escude, J. 2004, ApJ , 602, 1 

Ciardi, B., & Ferrara, A. 2005, Space Science Reviews, 116, 625 

Ciardi, B., Scannapieco, E., Stoehr, F, Ferrara, A., Iliev, I. T., & Shapiro, 

R R. 2005, MNRAS , in press ( astro-ph/05 1 1 623 1 
Dijkstra, M., Haiman, Z., Rees, M. J., & Weinberg, D. H. 2004, ApJ , 601, 

666 
Fan, X., et al. 2001, AJ , 122, 2833 
— . 2005, in press at AJ 1 astro- ph/05 120801 
Field, G. B. 1958, Proc. I.R.E., 467240 
— . 1959a, ApJ , 129, 525 

— . 1959b, ApJ , 129, 536 

Furlanetto, S. R. 2006, submitted to MNRAS tastro-ph/0604040t 

Furlanetto, S. R., & Loeb, A. 2002, ApJ , 579, 1 

—.2004, ApJ ,611,642 

— . 2005, ApJ , 634, 1 

Furlanetto, S. R., McQuinn, M., & Hemquist, L. 2006a, MNRAS , 365, 

115 
Furlanetto, S. R., & Oh, S. R 2005, MNRAS , 363, 1031 
— . 2006, submitted to ApJ 
Furlanetto, S. R., Oh, S. P, & Briggs, R H. 2006b, Physics Reports, in 

preparation 
Furlanetto, S. R., Zaldaniaga, M., & Hemquist, L. 2004, ApJ , 613, 1 
Gilfanov, M., Grimm, H.-J., & Sunyaev, R. 2004, MNRAS , 347, L57 
Gnedin, N. Y., & Hui, L. 1998, MNRAS , 296, 44 
Grimm, H.-J., Gilfanov, M., & Sunyaev, R. 2003, MNRAS , 339, 793 
Haiman, Z., Quataert, E., & Bower, G. C. 2004, ApJ , 612, 698 
Hirata, C. M. 2005, submitted to MNRAS ( astro-ph/0 507102t 
Hogan, C. J., & Rees, M. J. 1979, MNRAS , 188, 791 
Iliev, I. T., Mellema, G., Pen, U.-L., Merz, H., Shapiro, P R., & Alvarez, 

M. A. 2005a, submitted to MNRAS I astro-ph/05 121 87 1 
Iliev, 1. T., Scannapieco, E., Martel, H., & Shapiro, R R. 2003, MNRAS , 

341, 81 
Iliev, I. T., Scannapieco, E., & Shapiro, R R. 2005b, ApJ , 624, 491 
Iliev, 1. T., Shapiro, P R., Fen-ara, A., & Martel, H. 2002, ApJ , 572, L123 
loka, K., & Meszaros, R 2005, ApJ , 619, 684 
Kitayama, T., & Ikeuchi, S. 2000, ApJ , 529, 615 
Kuhlen, M., Madau, P, & Montgomeiy, R. 2006, ApJ , 637, LI 
Liszt, H. 2001. A&A , 371, 698 

Madau, R, Meiksin, A., & Rees, M. J. 1997, ApJ , 475, 429 
Malhotra, S., & Rhoads, J. E. 2004, ApJ , 617, L5 

Navarro, J. R, Rrenk, C. S., & White, S. D. M. 1997, ApJ , 490, 493 (NFW) 
Oh, S. P, & Haiman, Z. 2003, MNRAS , 346, 456 
Page, L., et al. 2006, submitted to ApJ I astro-ph/060345Ul 
Press, W. H., & Schechter, P 1974, ApJ , 187, 425 
Pritchard, J. R., & Furianetto, S. R. 2005, MNRAS , in press 

(astt-o-ph/0508381 1 
Ranalli, R, Comastri, A., & Setti, G. 2003, A&A , 399, 39 
Schaye, J. 2001, ApJ , 559, 507 
Scott, D., & Rees, M. J. 1990, MNRAS , 247, 510 
Seager, S., Sasselov, D. D., & Scott, D. 1999, ApJ , 523, LI 
Shapiro, P R.. Iliev, I. T., & Raga, A. C. 2004, MNRAS , 348, 753 
Sheth, R. K., & Tormen, G. 1999. MNRAS , 308. 1 19 
Smith, R J. 1966, Plan. Space Sci., 14, 929 



© 0000 RAS, MNRAS 000, 000-000 



The 21 cm forest 9 



Spergel, D. N., et al. 2006, submitted to ApJ (astro-ph/0603449t 

Taniguchi, Y, et al. 2005, PASJ , 57, 165 

Thoul, A. A., & Weinberg, D. H. 1996, ApJ , 465, 608 

Wouthuysen, S. A. 1952, AJ , 57, 31 

Zahn, O., et al. 2006, in preparation 

Zaldarriaga, M., Furlanetto, S. R., & Hemquist, L. 2004, ApJ , 608, 622 

Zygelman, B. 2005, ApJ , 622, 1356 



© 0000 RAS, MNRAS 000, 000-000 



